Voltage controlled exchange energies of a two electron silicon double quantum dot 
with and without charge defects in the dielectric 
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Quantum dots are artificial atoms used for a multitude of purposes. Charge defects are commonly 
present and can significantly perturb the designed energy spectrum and purpose of the dots. Voltage 
controlled exchange energy in silicon double quantum dots (DQD) represents a system that is very 
sensitive to charge position and is of interest for quantum computing. We calculate the energy 
spectrum of the silicon double quantum dot system using a full configuration interaction that uses 
tight binding single particle wavef unctions. This approach allows us to analyze atomic scale charge 
perturbations of the DQD while accounting for the details of the complex momentum space physics 
of silicon (i.e., valley and valley-orbit physics). We analyze how the energy levels and exchange 
curves for a DQD are affected by nearby charge defects at various positions relative to the dot, 
which are consistent with defects expected in the metal-oxide-semiconductor system. 
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Electrostatically defined quantum dots (QDs) in semi- 
conductors are one of the most promising systems for re- 
alizing a scalable quantum computer These systems 
provide the opportunity to engineer and control quantum 
mechanical properties through conventional electronics 
Q, and all elements of a qubit have been demonstrated 
in a GaAs double quantum dot (DQD) [3|. The presence 
of nuclear spins in GaAs, however, limit its spin coher- 
ence times to order of nanoseconds d, Q without any 
quantum control. Silicon, on the other hand, can have 
exceptionally long spin coherence times, often on the or- 
der of milliseconds to seconds 0, [H[ , and is therefore a 
very promising host material for solid-state qubits. 

After a decade of extensive research, experiments are 
beginning to build and measure few electron silicon QDs 
repeatedly [sHlH . The recent demonstration of single- 
shot spin readout in silicon also marks a significant mile- 
stone towards building a functional silicon qubit [lij . 
Several groups are now achieving coherent control ma- 
nipulation of silicon spins. However, charge defects rep- 
resent a potential challenge to forming predictable QD 
energetics and spin behavior. Inevitable charge defects 
in metal-oxide-semiconductor (MOS), for example, could 
strongly localize the electron even within a dot. In this 
letter, we investigate the effect of fixed charges in the 
dielectric on the voltage tuned exchange (J) curve of a 
silicon MOS DQD. We investigate the impact of strong 
localizing Coulomb centers on the J-curve tunability ac- 
counting for the full silicon band structure. 

We have developed an atomistic tight-binding (TB) 
based full configuration interaction (FCI) method to 
compute the multi-electron states of a DQD [15]. The 
single particle wavefunction is solved using the TB 
method for a DQD potential superimposed on the crys- 
tal potential. The QD wavefunctions from the TB solu- 
tion form the single electron basis for the FCI. We com- 
bine the atomic scale effects with an exact many electron 



method to capture excited electron configurations, ex- 
change, and correlation effects with unprecedented accu- 
racy. The method enables simulation of realistic devices 
as millions of atoms can be simulated in high performance 
computing clusters [16|, and in general, the atomistic TB 
method captures realistic details of the devices including 
miscuts, step roughness, alloy disorder, valley splitting, 
confinement geometries, strain, and applied fields. 
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FIG. 1: Binding energies due to (a) a single charge defect 
as a function of depth D into the oxide, and (b) two charge 
defects as a function of their separation R for two different 
depths D. Inset of (a) shows the probability of finding two 
defects a distance R apart for various defect densities. 

Fixed positive charges in the oxide near the Si-Si02 in- 
terface are common in MOS devices. The detailed chem- 
ical basis for the formation of these charge defects have 
been explored in Refs I?], [H| . Typical densities of these 
defects can range between 10 9 /cm 2 to 10 12 /cm 2 [l9| and 
positive charge is a common polarity of the fixed defects. 
Each defect charge produces a Coulomb potential, the 
tail of which penetrates into silicon and forms a shallow 
potential well that can bind at least one electron below 
the conduction band (CB). 

Fig. 1(a) shows the binding energy of an isolated defect 
in Si02 as a function of the defect depth D, schemati- 
cally shown in the top inset of Fig. 1(a), computed from 
the TB method. A single defect can bind an electron 
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in silicon with energies of a few meV, comparable to the 
orbital energies of a QD. Inset of Fig. 1(a) shows the 
distribution of defect separation distances for various de- 
fect densities, as given by a Poisson distribution, whereas 
Fig. 1(b) shows the binding energies of two defects as a 
function of their lateral separation R (inset of Fig. 1(b)) 
at a depth D into the oxide. This shows that a cluster 
of defects can bind electrons even more strongly and are 
detrimental to QDs. The binding energies compare well 
with recent measurements of defect states using electron 
spin resonance techniques [20|. For these calculations, we 
used a large simulation domain of 100 nm x 100 nm x 
50 nm to avoid finite size effects. 

The single electron TB Hamiltonian hi of a defect and 
a DQD, corresponding to the i th electron, is expressed as 



hi 



Hq + min{a(x — L) 2 + e, a(x + L) 2 } + ay 2 + eF z , 
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where Hq is the TB Hamiltonian of the host silicon 
formulated semi-empirically with the 10 band sp 3 d 5 s* 
model [Hj . The second and the third terms describe the 
2D parabolic potential energy of the QDs with curva- 
ture a and center-to-center separation of 2L. The left 
dot is subjected to a detuning of e relative to the right 
dot. The 4th term is the potential energy due to a uni- 
form vertical electric field F z that confines the electrons 
at the interface. The 5th term is the Coulomb poten- 
tial energy of a defect located at (xd,yd,—D) (interface 
being at z = 0) along with its image charge correction 
given by the factor Q = eSl e ^.° 2 21] . We have used 

(: Si | COlv_/2 

esi = 11.9eo and esio 2 = 3.9eo as the dielectric constants 
of Si and Si02 respectively, with eo as the permittivity 
of free space. Eq (pQ) is defined in the silicon region only 
with z > 0. The Si-Si02 interface is modeled as a hy- 
drogen passivated surface [22] . as used in earlier works 
[23| . The full Hamiltonian is solved with a parallel Block 
Lanczos algorithm to extract the relevant eigenenergies 
and wavefunctions near the CB minimum using the Na- 
noelectronic Modeling Tool (NEMO 3D) [l6|. The DQD 
system simulated in this work is comprised of about 1 
million atoms, and was solved typically on 40 processors 
in 10 hours. 

Using a set of lowest energy single electron states ob- 
tained from the TB Hamiltonian, all possible antisymm- 
teric two-electron configurations are constructed, and the 
two-electron Hamiltonian shown below is diagonalized in 
this basis, 



H = hi + h 2 
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(2) 



where h\ and h 2 are given by eq[TJ and the 3rd term is 
the electron-electron repulsion term with electron coor- 
dinates 7*1 and r 2 . The solution of eq (|2]) yields the two- 
electron states of the DQD-defect system. On average, 



we found about 12 single particle states corresponding to 
66 two-electron configurations were needed for conver- 
gence. 
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FIG. 2: (a) Low lying energy states (total energies) of a 
double quantum dot without any charge defects as a function 
of the detuning energy e indicate regions of e where the DQD 
is more (1,1) or (0,2) in nature. The dot has curvature a = 
0.0001 eV/nm 2 equivalent to a ID harmonic well with energy 
£0 = 8 meV. The vertical field is F z = 5 MV/m. The energies 
are with respect to the CB minima of bulk silicon at F z — 0. 
Inset: ID schematic of the detuned DQD potential, (b) Ideal 
(no defect) exchange (J) curves (i.e. various singlet-triplet 
splittings) of the DQD extracted from (a) as a function of e. 

In Fig. 2(a), we show the low energy two-electron spec- 
trum (i.e. total energy of the two-electron Hamiltonian 
of eq (2)) of the silicon DQD as a function of the detuning 
energy e (shown in the inset). At e = 0, the DQD is in 
a (1,1) charge configuration with each electron in a sep- 
arate dot. As e is increased adiabatically, the DQD goes 
through a gradual charge transition to the (0,2) state. 
Since the wavefunction symmetries and the spatial extent 
change during this bias sweep, the exchange (J) energy 
varies as a function of e. It has been demonstrated ex- 
perimentally in Ref Q for a GaAs DQD that the voltage 
controlled exchange can be used to provide a coherent 
rotation of the qubit encoded in the singlet-triplet basis. 

The bias dependence of the 2e spectrum is well- 
established for a GaAs DQD @, @|. In Fig. 2, we have 
extended this calculation to a silicon DQD, where the 
multiple conduction band (CB) valleys add an additional 
degree of freedom to the electronic states. The TB calcu- 
lations take into account all six CB valleys by represent- 
ing the full bandstructure of silicon. Due to the applied 
vertical E-field of 5 MV/m in the simulations, the two k z 
valleys are split from the other four CB valleys. Interface 
and E-field induced inter- valley coupling causes a further 
valley splitting of about 0.1 meV between the two lowest 
k z valleys for the passivated interface used in this work. 

The salient feature of the two-electron energies, shown 
in Fig. 2(a), is the presence of a multiplet of levels which 
vary with the bias in a similar manner, whereas as a 
single level would be expected in a GaAs DQD. These 
multiplets are a result of the two k z valleys, as there ex- 
ist multiple states with the same orbital envelope and 
spin wavefunctions, but with different rapidly oscillating 
Bloch wave components i.e. valley configurations. In the 
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case of Fig. 2(a), each multiplet consists of three distinct 
lines offset by a roughly constant valley splitting. The ex- 
change and Coulomb integrals used to evaluate the FCI 
Hamiltonian are strongly dependent on the valley config- 
uration of the states [24[, and methods that use a single 
valley approximation for silicon DQDs ignore potentially 
important details of the spin and valley physics. Since 
Fig. 2(a) is for B=0 T, the polarized triplets (T+ and 
T_) and the unpolarized triplet To are degenerate, and 
only the To triplets are shown. 

The exchange energy (J), defined as the splitting be- 
tween the lowest singlet and triplet states, can assume 
multiple definitions here due to the valley degrees of free- 
dom. In Fig. 2(a), we have labeled the energy difference 
of the lowest four triplets relative to ground state singlet, 
as JO, Jl, J2 and J3 in order of increasing magnitude. 
Fig. 2(a) shows that there is a low lying triplet almost 
degenerate with the first excited singlet. These are states 
with predominantly orthogonal valley character 24], 25| 
and hence small exchange splitting JO. The magnitude 
of JO primarily depends on the valley splitting caused by 
the interface and the vertical E-field. The Jl splitting 
is due to the lowest triplet with similar valley configura- 
tion as the ground state singlet, and is analogous to the 
J-curve of a GaAs DQD. The magnitude of Jl mostly 
depends on the curvature of the dots, the magnitude of 
the tunnel barrier, and their separation distances. For 
large valley splitting, the JO curve can be higher in en- 
ergy than the Jl curve. Recent experiments have shown 
that the relaxation between triplets of different valley 
configuration can be strongly suppressed [26|, and hence 
our primary focus in this work will be the Jl splitting, 
although we will also show how a defect perturbs the JO 
splitting. 

The four J-curves are plotted as a function of e in 
Fig. 2(b). The JO curve remains insensitive to detuning 
as the valley splitting is not affected by the lateral E- 
field. In the (1,1) charge configuration realized at small 
e, the high tunnel barrier between the dots reduces the 
overlap between the electronic wavefunctions. Hence the 
exchange energy Jl is small. As e is increased the (1,1) 
and the (0,2) singlets anti-cross, and the later evolves into 
the ground state. The triplet, however, remains in the 
(1,1) configuration as the (0,2) triplet is at higher energy. 
In this regime, Jl increases steeply with e as the singlet- 
triplet splitting increases (shown in Fig. 2(a) by the red 
arrow). At large enough e the (1,1) and (0,2) triplets 
eventually anti-cross each other, and the lowest singlets 
and triplets are all in the (0,2) configuration. Since the 
(0,2) configuration is a high overlap system with the 2e 
confined in a single dot, the exchange energy is high. Any 
further detuning of the dots has only negligible effect on 
the (0,2) configuration, and Jl becomes insensitive to e. 
J2 and J3 behave in a similar manner to Jl. 

Having established the voltage dependence of the J- 
curve in an ideal silicon DQD, we now analyze how the 
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FIG. 3: (a) Schematic of a DQD labeled L and R with oxide 
charges (i) at the tunnel barrier, (ii) at the center of one dot, 
and (iii) at various locations relative to the center of one dot. 
The corresponding exchange curves (b) Jl and (c) JO of the 
DQD for a charge defect in the oxide at the tunnel barrier 
between the dots (case (a)(i)). (d) Jl and (e) JO curves for 
an oxide defect at the center of the left dot (case(a)(ii)). (f) 
Jl and (g) JO curves for different defect locations relative to 
the center of the dot (case (a) (iii)). 



J-curve is perturbed in the presence of nearby oxide 
charges. In Fig. 3(b) and 3(c), we show the effect of a 
defect in the tunnel barrier between the dots (schematic 
of Fig. 3(a) (i)) on Jl and JO respectively. A positive 
charge defect in the oxide lowers the potential barrier 
between the dots, and increases the overlap between the 
electronic wavefunctions in the (1,1) charge configura- 
tion. As a result, Jl increases at low detuning energy 
relative to the ideal Jl-curve shown by the red dotted 
line. A lower barrier also makes the charge transition 
from (1,1) to (0,2) smoother, and the Jl-curve flattens 
out in the transition region. As shown in Fig. 3(b), in- 
creasing the charge magnitude or decreasing the defect 
depth impacts the potential barrier more, as the corre- 
sponding Jl-curves shift up and flatten out. Fig. 3(c) 
shows that the JO curve typically shifts up if the influence 
of the defect increases. This is because a stronger defect 



potential enhances the vertical E-field, which increases 
the valley splitting in both dots. 

As the positive defect charge is increased, the DQD is 
essentially transformed into a SQD insensitive to e. This 
means that if there are sufficient number of defects in 
close proximity to the tunnel barrier, the tunability of the 
Jl-curve is hampered. In such a case, attempts to form 
a DQD could be futile as the electrons will always reside 
in a SQD charge configuration. Another consequence of 
a defect in the tunnel barrier is that it could result in an 
always 'on' exchange gate, as the exchange interaction in 
the (1,1) occupation of a DQD is increased. 

In Fig. 3(d) and 3(e), we show the J-curves for a defect 
located at the center of the left dot in the xy-plane with 
various depths into the oxide (schematic of Fig. 3(a)(ii)). 
In this case, the defect lowers the potential of the left dot 
relative to the right even at e = 0, producing a natural 
detuning bias for the dots. The magnitude of this de- 
tuning increases as the defect is located closer to the dot 
(i.e. depth D decreases). Fig. 3(d) shows that the Jl- 
curves are translated more in e as D decreases, as more 
detuning bias needs to be applied to compensate for the 
defect induced detuning and to restore the electrons to 
the (1,1) charge configuration. Since the defect has only 
a small effect on the tunnel barrier, the slope of the J- 
curves remain the same, unlike the cases considered in 
Fig. 3(b). The presence of a defect in one dot relative 
to the other also causes an asymmetry in the J-curve 
between the (0,2) and the (2,0) configurations as the de- 
fect may produce a stronger confinement for the electrons 
and may result in a stronger exchange splitting. The J0- 
curves of Fig. 3(e) exhibit a slope with e due to the fact 
that the defect produces a larger valley-splitting in the 
left dot compared to the right. Hence the JO curve de- 
creases in magnitude as the electrons go from a (2,0) to 
a (0,2) transition. 

In Fig. 3(f) and 3(g), we show the effect of a defect 
at various locations relative to the dot (schematic of Fig. 
3 (a) (hi)). The defect depth D is fixed at 2 nm, and the 
dot centers are 20nm apart along the x-axis. The label 
x = +5 nm indicates a defect in between the left dot and 
the tunnel barrier, displaced 5 nm from the left dot L 
in the x-direction. Similarly, the label x = — 5 nm in- 
dicates a defect located 5 nm right of the center of the 
dot and 15 nm from the tunnel barrier. The J-curves 
corresponding to these two cases can be understood as 
a combination of Fig. 3(b)-(e). The defect essentially 
detunes the left dot and thus translates the Jl-curve, 
but to a lesser extent than it does in Fig. 3(c) where 
the defect is at the dot center. A defect in between the 
dot and the tunnel barrier, in addition, lowers the tunnel 
barrier, which manifests as a flatter Jl-curve. The label 
y = +5 indicates a defect displaced 5 nm from the dot 
center perpendicular to the DQD axis. The Jl-curve in 
this case is also translated, and lies between the two pre- 
vious curves, as expected. Similarly, the JO-curves also 



exhibit an increased magnitude and a slope as expected. 

We have done a full configuration interaction calcula- 
tion of a Si DQD using TB wavefunctions that account 
for valleys and atomistic details. Our results show that 
charge defects in the oxide near a silicon DQD can de- 
termine the characteristics of the multi- valley exchange- 
curve. Defects in the tunnel barrier affect the J-curve 
most as they can hamper the tunability of a SQD to a 
DQD. Defects at other locations mainly manifest as a 
translation in the J-curve to a different detuning bias. 
The combined TB and CI method developed to perform 
these calculations represent a significant advancement 
in the computational simulation methods for nanostruc- 
tures, as it integrates the atomic scale details with an 
exact many-electron theory, and can easily scale up to 
million atom systems. 
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